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Abstract 


Quantitative criteria are presented for model simplification, or order reduction, such that the reduced order 
model may be used to synthesize and evaluate a control law, and the stability and stability robustness obtained using 
the reduced-order model will be preserved when controlling the full-order system. The error introduced due to model 
simplification is treated as modeling uncertainty, and some of the results from multivariable robustness theory are 
brought to bear on the model simplification problem. A numerical procedure developed previously is shown to lead 
to results that meet the necessary criteria. The procedure is applied to reduce the model of a flexible aircraft. Also, 
the importance of the control law itself, in meeting the modeling criteria, is underscored. An example is included 
that demonstrates that an apparently robust control law actually amplifies modest modeling errors in the critical 
frequency region, and leads to undesirable results. The cause of this problem is identified to be associated with the 
canceling of lightly-damped transmission zeroes in the plant. 

Whether the engineer is developing a system model for dynamic analysis, control law synthesis, or 
simulation, a simple low-order model with the requisite validity is desirable for a variety of practice reasons. The 
question arises, therefore, as to how to obtain such a simple yet valid model. Even more fundamental is the question 
of what model characteristic are important such that one may strive to retain them. Although the initial question 
has been addressed for some time, from the attention still paid to model and controller order reduction (c.f. Refs. 1,2), 
it appears that the issues still remain unresolved. 

In Refs 3-6, some previous offerings on the subject are presented. In this paper, discussion will continue, 
in the attempt to expand on some of the earlier results, to further clarify the theoretical basis behind the proposed 
methodology, and to reveal some important aspects of not only model-simplification, but also control-law synthesis 
for elastic vehicles. 


1 , Criteria for Modeling 

The objective in model simplification, as with all system modeling, is to develop a fundamental 
understanding of the system in question. For the model to be useful, it should predict to the required engineering 
accuracy the behavior of the actual system. Note that it does not have to predict with perfect accuracy, and that is 
not possible anyway. The required accuracy depends on the application for which the model is intended. 

In this paper, as in Refs. 3 - 6, the intended application of the model is to predict the behavior of the 
system when it is subject to feedback action, as shown, for example, in Fig. 1. Clearly, then, the critical 
characteristics of the actual system that must be adequately captured by the model are those characteristics important 
in a feedback system. (Note that the feedback action could represent an automatic control system, as well as that of a 
human, or manual controller.) Finally, the existence of a sufficiently valid, although perhaps complex model for the 
system is assumed to be available - admittedly a big assumption. Further, if this model is infinite-dimensional 
and/or non-linear, it is assumed that a locally linearized, finite-dimensional model may be obtained. The original 
(complex) model will be denoted as Q, while the linear model will be denoted as G. 

As a result of any simplification process, differences between the more-accurate model and the simple model 
arise. Or conceptually, if Gr is a simpler model for G, the model-simplification error may be considered to be AG = 
G - Gr. These errors are key to the research presented here. In contrast, model-simplification errors arising due to 
the development of G, or A Q = Q - G , will be considered only indirectly. 


650 


The critical question then is what errors AG are critical, or should be minimized, and what procedure will do so? 
The answer to the first part of the question could be that AG’s critical in a feedback loop should be minimized. Further, 
if these AG's are interpreted more generally as model uncertainty rather that model-reduction error, the recent research on 
multi-variable robustness theory may be bought to bear on the model-simplification problem. This is the main idea in 
this research. 


2. Robustness and Model Reduction 

In this section, some key results from robust control theory will be noted, and they will be interpreted in the 
context of the model reduction problem. 

With reference to the system shown in Fig. 2, G R is the transfer-function-matrix representation of this 
simplified model, AG(s) in the analogous representation of the model-simplification error, and the full-order linear model 
is G = G R + AG. Likewise, K(s) is the matrix of control compensators, perhaps to be designed using G R . Clearly in 
this context, one desires that the K(s) so obtained will control the "true” G(s) as predicted through the use of G R . 
Attention is now turned to exposing the critical AG's via multi-variable Nyquist theory.™ 

Let d>(s) be an analytic function of the complex variable s, and let the number of zeros of d>(s) in the open right 
half of the complex plane be denoted as z. Then the Principle of the Argument states that 

N 

R — >oo (O, d*(s), D r ) — z 

or the number (N) of clockwise encirclements of the origin made by the image of the contour D R , under the mapping of 
<f>(s), as s travels clockwise around D R , equals z. Here D R is the "Nyquist D contour" that encloses the entire right-half 
of the complex plane. Clearly, with regards to stability, the <f>(s) of interest is the closed-loop characteristic polynomial 
of the feedback system, denoted by O cl (s). 

Now, as shown in Ref. 8, and elsewhere, and referring to Fig. 1, for example, 

dVCs) = <t> OL (s) det [I + GK] 

= <b OL (s) det [I + KG] 

where <f> OL (s) is the characteristic polynomial of the open-loop system KG(s) or GK(s) 
function matrix GK(s) or KG(s) has the state-space realization 

X = ^GK X+ ®GK e 

y = CgK X 

then d> OL (s) = det[sl - A CIC ], and the zeros of <J> 0L (s) are the open-loop poles of the system. Note that Eqn. 1 may therefore 
be re-written as 


( 1 ) 

. That is, if either the transfer 


<X> OL (s) = det[sl - A gk ] det [I + C GIC [si - A GK ]'' B GK ] 


Now if the number of right-half-plane zeros of <t> 0L (s) is p, then the number of right-half-plane zeros of det [I + 
C GK [si - A 01C ]'’ B otc ] must be z-p. Furthermore, from the Principle of the Argument 

N 

R->oo (O, det [I + C GK (si - A gk ) 1 B 01( ) ,D r ) = z-p 


Consequently, if p is known, z may be deduced from 


z = p + (z - p) 

= p + [r— joo (O, det [I + C 0K (si - A^) 1 B ok ] ,D r )] 
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or closed-loop stability is determined from knowledge of p and the examination of the Nyquist contour for det[I + GK] or 
det [I + KG]. Therefore, the closed-loop system is stable if and only if the Nyquist contour for det[I + GK](= det[I +KG]) 
encircles the origin counterclockwise exactly p times. 

Of course the determination of z is possible from other means, and the real utility of the above fact is in 
defining the concept of relative stability, and in identifying factors that are critical to closed-loop system stability. These 
issues are of special import here. 

Consider the model error, or uncertainty, to be AG (as in Fig. 2), and assume that K is such that KG R leads to a 
stable closed-loop system with good stability margins. (Note this assumption should always be true as it involves a key 
objective in determining K(s) using G R to begin with.) Then if (assumption 1) the number of right-half-plane poles of 
KG (= p) is identical to the number of right-half-plane poles of KG R (= Pr), K will stabilize G if and only if (assumption 
2 ) 


R^oo (0, del [I + GK], D R ) = R ^oo (O, det [I + G R K], D R ) 

or the number of encirclements of the origin made by the Nyquist contours associated with G and with G R are identical. 

Stability is guaranteed as follows: 

Let z = no. of unstable closed-loop poles of the KG loop. 

Zr = no. of unstable closed-loop poles of the KG R loop 
p, p R = defined above 

Then to show stability (or z = 0), note that if (assumption 1) p R = p, then 


Z = Z R - (2* - Pr) + ( Z - P) 

By the assumption KG R leads to a stable system, Zr = 0, ancl from assumption 2, (z R - p R ) = (z-p). Hence, z = 0. 


This now establishes in a meaningful way, qualitative criteria for model simplification, the simplification must 
at least lead to AG’s such that assumption 1 and 2 are satisfied. But the criteria goes further. Not only must stability of 
the KG loop be assured (i.e., z = 0) but the margins "designed" into KG R should carry over to the closed-loop system 
associated with KG. Otherwise, the K so designed would not be satisfactory. It is for this reason that any model 
reduction technique that just assures stability of the full-order closed-loop system may not be good enough! 

To satisfy assumption 2, or to assure that the number of encirclements of the origin is unchanged due to AG, 
requires that {9] 


det [ I + G r K + eAGK] * 0 v co>0, e e [0,1] (2) 

In other words, if as the Nyquist contour for det[I + G R K] is continually warped to that for det [I + GK] the origin is 
never intersected, the number of encirclements of the origin cannot change. Furthermore, Eqn. 2 is assured if (c.f.. Ref. 
9) 

ft (AGK) <g[l + G r K] vg)>0 (3) 

Finally, it is known that an alternative to Eqn. 3 is 

0%) <<?[! + (GrK) 1 ] = g {[G r K (I + GrK) 1 ] 1 } v co>0 (4) 

where E m = Gr^AG 
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The above expressions (Eqns. 2 - 4) may be extended by breaking the frequency domain (0 < co < oo) into the 
domains (0 < co < co*) and (co* < co < oo). Note that these domains are non-intersecting. Now it can be argued that 
Eqn. 2 will be satisfied if 

det [ I + G r K + eAGK] * 0 (0< to < to*) (5) 

(0 < e < 1) 

and 

det [ I + G r K + eAGK] *■ 0 (to* < to < oo) (6) 

(0 < e < 1) 

Further, Eqn. 5 is assured if Eqn. 3 is satisfied for co < co*, while satisfying Eqn. 4 for co > to* assures that Eqn. 6 is 
satisfied. Hence, in such a situation, Eqn. 2 is satisfied. 

By Eqns. 3 and 4, quantitative criteria on critical AG's are established. Further, the overall strategy for model 
simplification becomes apparent, and the interaction between model simplification and control law synthesis is 
underscored. Regarding the later, it should be clear that the allowable AG's (those that do not destroy closed-loop 
stability of the full-order system controlled by K(s)) depend on K itself. In other words, designing a "good" K(s) 
increases that allowable AG, while designing a bad one may put very strict limitations on the allowable AG, and hence 
model accuracy. The former K(s) is robust, the latter is not. 

Regarding the model simplification strategy, then, first observe the right side of Eqn. 3. When ct (G r K) »1, 

ct [I + G r K] = ct (G r K). Conversely, when ct(G r K)«1, ct [I + G„K] = 1. Finally, the CT [I + G R K) will take on its 
minimum value in the frequency range where ct, (G r K) = 1. The frequency range where the latter occurs is of course the 
(multi-variable) gain crossover region. Consequently, it is this frequency range where the AG must be the smallest, and 
this can be assured if each element of the AG matrix is small in this frequency range. 

Also, noting the above discussion, Eqn. 3 may be satisfied by rather large AG in any frequency range where ct [I 
+ G r K] is large, and this will occur when ct (G r K) is large. If K is designed to give a good classical Bode loop shape, ct 
(G r K) will be large for frequencies below crossover. 151 

Now consider Eqn. 4. When ct (G h K) «1, ct (G r K)''» 1, and ct {I + (G R Ky ! ] = ct (G r K) - 1 »1. Hence the 

allowable AG may also be rather large in this case. Further, if K yields a good loop shape, or is well attenuated, at high 
frequencies, a (G r K) will be small for frequencies above crossover. So clearly, the AG must be smallest in the region of 

multi-variable crossover, while if K yields a good bode loop shape, rather large AG elsewhere may be acceptable and 
Eqns. 3 and 4 may be satisfied. The above discussion is summarized in Fig. 3. 

The final issue to be addressed is that of satisfying assumption 1, or the number of unstable poles of KG R must 
be identical to the number of unstable poles of KG. First note that this is equivalent to requiring the number of unstable 
poles of G and G R to be the same, since only one K is involved. Then observe that the poles of G are the poles of G R + 
AG, which consists of the poles of G R plus the poles of AG. Hence to satisfy assumption 1, AG must be stable. 

Attention will now turn to some additional criteria arising from performance considerations rather that from 
stability robustness. The system to be considered is that shown in Fig. 4. The vector of responses Y(s) is given by 

Y = [I + (G, + AGj)K] -1 (G, + AG,)K(Y C - N) 

+ a + (G, + AG,)K]-> (G 2 + AGj) D 

Here G, is the reduced-order model for the response of G to control inputs, where G 2 is the reduced-order model for the 
response of G to disturbances being considered. AG, and AG 2 are the analogous model-simplification errors. 
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The first observation to be made is that stability and stability robustness depends on G l and AG P not on G 2 and 
AG 2 Note that the poles of (G l + AG t ) are the poles of the "true" plant G, as are the poles of G 2 + AG 2 . Hence if K 
stabilizes G, which will be assured if G { and AG { satisfy the criteria developed previously, K must therefore stabilize (G 2 
+AG 2 ). This is significant since some (stable) poles of G may be approximately cancelled by some zeroes for the 

transfer functions governing responses to control inputs, but not cancelled in those governing responses to disturbances. 
Cancelling these poles to obtain G 1 has raised questions by some as to whether those poles so cancelled could lead to 
problems later in analysis. The answer appears to be that they will not if G 2 is obtained such that those poles are 
retained. But from the above discussion on stability, the only reason to keep these poles in G 2 (that by assumption are 
not approximately cancelled) Is such that the disturbance -rejection performance predicted using G 2 (when designing K, for 
example) will be reasonably accurate. 


Finally, noting that the disturbance response due to AG 2 is 

Y Dj = [I + (Gj + AGjjK]* 1 AG 2 D 

for good performance prediction (Y D ^ small), AG 2 should tend to be small whenever D is large and (Gj + AGJK is small. 

But here again, if K is designed to obtain a "good loop,” it will be designed such that G X K (and by implication (G t + 
AGj) K) will be large over the frequency range where D is large. Consequently, this should not pose stringent 
requirements on AG 2 . 

In ending this section, it is worth noting that assuming K is designed properly has been critical. By doing so, 
one lakes advantage of one of the basic advantages of a good feedback system, reduction in sensitivity to plant (or plant 
model) variations. This allows the development of a modeling procedure that focuses on the really critical problem of 
obtaining a good model in the crossover region. 


3. Methodology and Sample Results 

The procedure offered was discussed in detail in Ref. 5, and the computational technique is summarized again in 
Table 1. The technique is a frequency weighted internally-balanced approach, with stable factorization in the case of an 
unstable plant G. The stable factorization procedure sets the unstable subsystem of G aside via partial fraction 
expansion, leaving the remaining subsystem G $ stable. This stable subsystem is then reduced, such that a stable reduced 
order model G RS is guaranteed. The unstable subsystem is then rejoined with G ^ to obtain the final reduced-order model 
G r . By this procedure, the number of unstable poles of G are preserved. In fact the unstable poles in G are exactly 
retained in G R . 

The internally balanced technique 1101 requires the frequency-weighting extension 151 since the basic technique leads 
to small model-simplification errors AG where the elements of G have large magnitude, which is not necessarily the 
crossover region. Further, a very poor model may be obtained where the elements of G have small magnitude. As will 
be shown later, this can be totally unacceptable. 


In Ref, 1 1, a frequency- weigh ted approach was also suggested, but the weighting required the knowledge of the 
compensator K, obtained using the full-order plant. Since designing a simple K using the simpler plant G R is the typical 

design objective, the above weighting is undesirable. In Ref. 5, it was noted that simply adding a weighting filter 
obtainable by inspection of the Bode plots of G and knowledge of the desired crossover frequency range let to excellent 
results. This filter is easily discarded after G R is retained. In the example presented later, it will be shown that this 
approach again appears quite acceptable. 

The key to the concept is the knowledge of the fact that the internally balanced approach yields a small AG 
where the elements of G have large magnitude. Heuristically, if a filter W(s) is used such that W(s)G(s) has large 
magnitude in the required frequency range, and if WG reduced such that WG R is obtained, then G R will have the desired 
properties. 
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Table 1 . Frequency Weighted Internally Balanced Reduction 


Given: System State space description A, B, C and weighting filter state space description A w , B w , C w . 

Find: r* order system 

Step 1: Solve for X and Y 


u bci r x x ] 2 7 r x x,n r a* oj + r 0 o] 

(P z w j J^X 21 X 12 y |X 21 X 22 J JC^B T A^J [o B w BjJ 

fpj o'] fr y 1 f" Y Y ,1 Pa Bcl fccr oj 

[pwB T AXJ [y 21 y 2 J [y 2] y 22 J [ 0 A w j [o 0 J 


Step 2: Find T and I where XY = TOT 1 , T = [T r T n _ r ], T T = [U r , U n _ r ] 


-o 


^ d 


r~ — i 

I 2 0 

I? = 

0 x n . r 


where 


Sr = diag(v ci v oi ) i = 1 r 

L n-r = diag( v ci v oi> * - r + 1 n 


v ■ v * ^ > v v > 0 

V CI V 01 ’ * Y C1 V 01 - v 


Step 3: r ^ order system is 


A r = UjATr 
B r = UjB 
C r = CT r 


As the example, consider an elastic aircraft identical to the configuration investigated in Refs. 3 and 6. This 
configuration is of reasonably conventional geometry with a low-aspect ratio swept wing, conventional tail, and canard. 
A numerical model for the longitudinal dynamics is available from the above references. Both rigid-body modes and four 
elastic modes (resulting in a ll lh order model) are included. The in-vacuo vibration frequencies are 6.3, 7.0, 10.6, and 
1 1.0 rad/s, and are representative for a supersoni^ypersonic cruise vehicle. These frequencies, furthermore, are all near 
the anticipated frequencies at crossover for the control systems to be designed. 

Control inputs are elevator deflection 5 E and canard deflection 5 C , while the disturbance is the perturbation in 
angle of attack due to atmospheric turbulence a g . Selected responses are vertical acceleration a z ’ measured at the cockpit 
and pitch rate q measured at the antinode of the first bending mode. Therefore, the flight and structural mode control 
loops in the context of Figure 4, might correspond to the following, for example 

Y = [a z f q] T 

U = [5 e 5 c F 

D = “g 

Obtaining the reduced order model G, was the subject of Ref. 6. An anticipated crossover frequency range (for 
G,K) was assumed as 1 to 10 rad/s. In that reference, it was also noted that a fourth-order for G, was sought based on 
the observation that the full order model has two oscillatory models in this frequency range. 
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Attention is now turned to the requirements for G r As a realistic example, the Dryden gust spectrum for 
turbulence is used to describe the disturbance. A fourth-order model for G 2 is sought based on the observation that the 
full order model has two oscillatory models in the frequency range where the spectrum of D is largest. This frequency 
range is coincidentally also 1 to 10 rad/s. 


The reduced order models for G x and G 2 were then obtained simultaneously from the frequency-weighted 
internally-balanced reduction technique [ 51 which was specifically developed to meet the criteria in Section 3. The 
frequency- weighting filter used was a band pass filter of unity magnitude in the 1 to 10 rad/s frequency range with 40 
db/dec roll off on either side of this frequency range. 


Table 2 contains the reduced order state space matrices A, B, C and D. Figures 5 through 10 show the reduced 
order and full order frequency response magnitudes for G x and G 2 . Observe that the reduced order model accuracy 
approximates the full order model in the 1 to 10 rad/s frequency range as desired. To complete this example, a simple 
control law, consisting of three constant gains was synthesized using the model G r The synthesis objective was to 
augment the damping of the first aeroeiastic mode with acceleration feedback to the canard, to augment the short period 
damping with pitch-rate feedback to the elevator, and to provide some response decoupling with a cross feed from the 
elevator to the canard. The resulting control law is of the following form 



Actuation effects were modeled with simple first-order lags, with comer frequencies at 20 r/s for both the canard and the 
elevator. 


Tab le 2 Order Model 


-.9932 

-2.013 

-5.593 

4^34 

.8294 

-.0137 

-.6638 

.2098 

-.0138 

.0121 

-.3175 

1732 

A IB 

c |d 

-.0507 

.0329 

-9.658 

-.5171 

-31.67 

35.92 

-593.7 

-281.4 

14.48 

-21.42 

-420.0 

-175.2 

13.59 

-24.38 

700.1 

342.5 

.0665 

-.03471 

.0017 

.0015 

0 

0 

0 

8.762 

.7218 

.9287 

-2.038 

52.01 

-244.5 

333.0 

y = q 

(r/s) 

u = 

A ( rad ) 


D = a (rad) 

az' 

(ft/s) 


5 C (rad) 





Shown in Fig. 1 1 is the plot of Eqn. 3, while Eqn. 4 is shown in Fig. 12. Note that although this control law 
did not result in high gain (large G, K) at low frequencies, Eqn. 3 was still satisfied below crossover region. Conversely, 
Eqn. 4 is satisfied, although barely, in the frequency range above crossover. Hence, from the argument in Section 2, if 
co* in Eqns. 5 and 6 is in the crossover region, stability is assured. For reference, the pitch-rate to elevator transfer 
function . .. - ^ ' 

q(s) _ 50 (0.33)[.13,4.84][.01,10.6][.03,11.0][.21,13.](45.) 

5 Ec (s) [.53,1.81][.15,4.78][.02 I 10.8][.03,11.][.19,13.3](19.)(69.) 


4, An Additional Criteria 

As noted in Section 3, the AG^ arising from the model simplification must satisfy stringent criteria Jn the 
crossover region, and if Eqn. 3 and/or 4 (or 5 and 6) is satisfied, closed-loop stability is assured. To be discussed here is 

the fact that the controller K should not be such that small AG is amplified such that a (AGK) becomes large. It will 
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be shown by example that this can easily occur where the magnitudes of G (or of the gy s) are small. Hence, the 
example will demonstrate why obtaining a good model in this situation is important (recall that unweighted balanced 
reduction has a problem here), and some implications regarding control-law synthesis will also arise. 

Consider the simple scalar plant 

g(s) = (s 2 + .04s + l 2 ) 

s(s 2 + .032s + 0.8 T 

The plant is stable and minimum phase, so a robust control law should be obtainable. Using LQG/LTR or H^, for 
example, the following compensator could be obtained. 

k(s) = 8(s 2 jr ,032s + 0.8 2 ) 

(s 2 * + .04s + l 2 ) (s + 8)' 

It can be easily verified that the loop shape kg is very good, yielding infinite gain margin, 90 degree phase margin, and 
good roll off above 8 rad/s. 

Now assume that the "true" plant is 

gtrue = 0.69 ( s * + • 048s + 122 > 
s(s 1 + ,032s + 0.8 7 ) 

or the numerator "frequency" is in error by 20% (1.0 -» 1.2). Note that this could occur, for example, if a vibration 

mode shape was slightly off in the modeling. Shown in Fig. 13 is the plot of Eqn. 3 for this example, and clearly a 
(Agk) > <r (1 + gk) at 1 r/s (the designed crossover frequency). Further, a quick check would show the kg^ loop to be 
unstable. But the lAgl = Ig - g true l (not shown) would be found to be rather modest at co = lr/s, with much larger lAgl at 
lower frequencies. The problem could be interpreted as one of the control law k amplifying the lAgl at co = 1 r/s, and this 
is confirmed from the plot of lk(jco)l in Fig. 14. 

Stability of the kg lfue loop would result, and Eqn. 3 satisfied, if the lk(jco)l at co = 1 rad/sec were simply 
reduced. This is accomplished with the following compensator 

k m0 d(s) = 8(s z + .032s + 0.8 2 ) 

(s 2 + 1.2s + l 2 ) (s + 8) 

or the damping of the complex compensator poles is increased, and the plant model zeroes close to the imaginary axis are 
not exactly cancelled. Clearly the loop shape with this compensator is not as "optimal" as the original, but this control 
law is more robust against this Ag. 

Noting that the problem arose with a modeling error that is associated with lightly-damped zeroes, the critical 
Ag was at a frequency (co = 1 r/s), where lg(jco)l was relatively small as shown in Fig. 15. Hence, obtaining a good 
model at this frequency is important. Furthermore, by attempting to cancel those lightly-damped zeroes in the plant, the 
original controller was very sensitive to their location. Increasing the damping of the compensator poles, as in a 
classical notch filter, made the loop more robust against the uncertainty in the location of these plant zeroes. 
(Incidentally, this can be accomplished with a modified LTR procedure, as noted in Ref. 12 and in another paper in 
preparation.) 

As a final remark, it is observed that lightly-damped zeroes in the compensator are different from similar zeroes 
in the plant since through the design and implementation of the compensator, the location of its zeroes may be more 
accurately defined. 


5. Conclusions 

Quantitative criteria are presented for model (or controller) simplification. The reduced order model (or 

controller) must well approximate the full-order system in the (multivariable) crossover region for stability, and stability 
robustness, to be assured. Bounds on the model-simplification error were noted, and if the bounds are satisfied, stability 
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is assured. It was also noted that the model reduction criteria were functions of the control law, and by synthesizing a 
robust control law, the criteria could be easier to satisfy. 

A numerical procedure, consisting of stable factorization with weighted balancing of coordinates has been 
shown, by example, to meet the above criteria. The example involved reducing an eleventh order linear model of an 
clastic aircraft to obtain a fourth-order model leading to the desired six transfer functions. 

Finally, another example demonstrated the importance of obtaining good agreement between the full- and 
reduced-order model in the crossover region, even where the transfer function (or functions) have relatively small 
magnitude. Furthermore, the^xample demonstrated that an apparently robust controller could in fact amplify small 
errors, and lead to unstable results. The problem would occur with any control law that had the effect of cancelling 
lightly-damped transmission zeroes of the plant model. 
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